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Theoretical  Aerodynamics  Note  ijfo.  33 


ANALYTICAL  INVESTIGATION  OF  FLOW  DURATION  IN  THE  CONVAIR  HYPERSONIC 
SHOCK  TUNNEL 


SUMMARY 


V 


>.7ne  analytical  procedure  used  to  calculate  the  steady  flow  duration  in  a 
hypersonic  shock  tunnel  is  presented i  The  results  of  this  procedure  applied  to 
the  Convair.  Hypersonic"  "Shock  Tunnel  are  given.  It  is  shown  that  for  very  low 


-v.  < 


diaphragm  pressure  ratios,  p! /pf'  .<,5^0,  the  steady  flow  duration  is  non 

•  7  J*- 

existent..  For  low  diaphragm  pressure  ratios,  ?40*,<  P^'/P^  ,<  2000,  which  in¬ 
cludes  the  tailored  interface  condition  the  steady  flow  duration  is  quite  small. 
This  is  a  result  of  pressure  attenuation  caused  by  the  expansion  waves  reflected 
from  the  end  of  the  driver.  I  f 
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SYMBOLS 


Velocity  of  sound 


Dimensionless  internal  energy  ratio 
Equivalence  factor  equation  (2) j 
Mach  number 


pressure 


Dimensionless  pressure  ratio 
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Particle  velocity 

\ 

Dimensionless  parti c?. 3  velocity  — 

SJ 

Position  along  shock  tube  measured  from  diaphragm 

c 

Ratio  of  specific  heats, 

v 

7  +  1 
7  -  1 

7-1 
2  7 

Density 

Rarefaction  wave 


Shock  wave 
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INTRODUCTION 

The  hypersonic  shock  tunnel  is  a  development  of  the  basic  shock  tube'''  in 
vhich  the  shock  tube  is  used  to  produce  a  reservoir  of  high  temperature,  high 
pressure  air  vhich  is  expanded  thru  a  nozzle  to  obtain  the  desired  test  condi¬ 
tions.  Even  though  the  shock  tube  is  basically  an  unsteady  flov  device  there 
are  several  regions  in  vhich  the  flov  is  steady  and  can  be  used  as  the  reser¬ 
voir.  The  durations  of  these  steady  flov  regions  in  the  basic  shock  tube  are 
quite  short,  being  of  the  order  of  a  millisecond  or  less. 

Tiie  steady  flov  duration  of  a  hypersonic  tunnel  is  determined  by  the 
duration  of  the  steady  state  reservoir.  The  required  flov  duration  may  be 
about  a  millisecond  for  adequate  measurement  of  pressure  and  heat  transfer 
characteristics  but  should  be  tvo  or  three  milliseconds  for  force  measurements. 

Two  modifications  of  the  basic  shock  tunnel  that  have  been  developed  to  increase 

2 

the  available  test  time  ere  the  reflected  shock  design  and  the  tailored  in- 
3 

ter face  design. 

The  duration  of  the  steady  flov  reservoir  and  hence  the  test  time  avail¬ 
able  is  a  function  of  facility  size  and  design  as  veil  as  of  the  characteristics 
of  the  gases  and  initial  pressure  ratio  used.  This  note  vill  explain  the 
analytical  procedure  used  in  calculating  steady  flow  duration  and  vill  present 
the  results  of  the  flov  duration  analyses  of  the  Convair  Hypersonic  Shock 
Tunnel . 
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Discussion 

Wave  Diagrams  and  Run  Time 

One  of  the  most  useful  devices  in  analyzing  shock  tube  flov  is  the  vave 
4  5 

diagram,  ’  The  wave  diagram  is  a  plot  of  the  wave  patterns  in  the  x,t,  plane 

where  in  the  case  of  a  shock  tube,  x  is  the  longitudinal  distance  along  the 

tube  and  t  is  the  time.  Figure  1  presents  the  wave  diagram  for  a  simple 

shock  tube  and  the  corresponding  flow  picture  inside  the  tube  at  a  time  t^. 

Figure  2  compares  wave  diagrams  for  the  simple  shock  tunnel,  the  reflected 

# 

shock  modification  and  the  tailored  interface  technique.  The  available  test 
time  for  each  of  these  methods  is  also  shown  in  Figure  2. 

The  simple  shock  tunnel, ^  merely  has  an  expansion  nozzle  at  the  down¬ 
stream  end  of  the  driven  section.  The  available  run  time  is  basically  the 
time  difference' between  the  arrival  of  the  initial  shock  wave  and  the  contact 
surface  at  the  nozzle. 

The  most  common  modification  to  the  basic  shock  tunnel  is  the  reflected 

shock  design  in  which  the  expansion  nozzle  is  replaced  by  a  conventional  con- 

2 

vergent- divergent  nozzle.  The  nozzle  acts  like  a  plate,  almost  completely 
reflecting  the  initial  shock  wave.  The  run  time  for  this  modification  is  the 
time  difference  between  the  arrival  of  the  initial  shock  vave  and  the  wave  re¬ 
flected  from  the  contact  surface. 

The  tailored  interface  shock  tunnel  is  merely  a  special  case  of  the  re- 
fxected  shock  modification. J  In  tills  case  the  acoustic  impedance  is  matched 
at  the  intersection  :f  the  reflected  shock  and  the  interface  so  that  no  wave 
is  reflected  from  the  interface.  The  available  test  time  is  the  time  difference 
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between  the  arrival  of  the  initial  shock  wave  and  the  unsteady  expansion  waves 
propagating  downstream  from  the  diaphragm. 

These  definitions  of  available  test  time  assume  that  the  expansion  waves 
propagating  upstream  and  reflecting  from  the  end  of  the  driver  do  not  overtake 
the  waves  propagation  downstream  before  the  events  described  above  take  place.' 

Wave  Diagram  Construction 
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The  method  for  constructing  a  wave  diagram  will  be  described  for  a  shock 
tunnel  with  an  area  change  at  the  diaphragm  and  a  convergent- divergent  nozzle 
at  the  end  of  the  driven  section.  The  wave  diagram  of  such  a  tunnel  is  shown 
in  Figure  3*  Also  shown  in  Figure  3  &re  the  pressure  distribution  and  wave 
picture  along  the  tube  at  a  time  t. 

The  first  stepA in  constructing  the  wave  diagram  is  a  knowledge  of  the 
initial  conditions  in  the  driver  aiw.  driven  sections  of  the  shock  tube  before 
rupture  of  the  diaphragm.  These  include  pressure,  temperature  speed  of  sound 
and  the  ratio  of  specific  heats,  7.  Given  these  initial  conditions  the  pres¬ 
sure  rise  across  the  initial  shock  can  be  determined  by  an  iterative  procedure 
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Since  A,  is  the  minimum  section,  I<L  is  equal  to  unity  but  M_  must 

^  ^a 

still  be  determined  in  order  to  solve  for  g  in  equation  (2) .  M~  can  be  deter- 

\ 

mined  from  the  folloving  equation  since  the  area  ratio  -7-  is  known. 

74+1  \ 


A 

\  _^b 


2  +  (yh-i)u  ‘ 

_ fia 

L2  + 


(3) 


For  the  Convair  facility  these  equations  have  been  solved  for  two  initial 
driver  conditions  and  the  variation  of  initial  shock  pressure  rise  as  a  function 
of  diaphragm  pressure  ratio  is  presented  in  Reference  8.  This  curve  is  repro¬ 
duced  herein  as  Figure  4. 

P2 

Given  the  pressure  rise  —  across  the  shock  the  other  quantities  across 

P1 

the  shock  can  be  obtained  from  the  cl  .arts  of  Reference  9.  The  two  quantities 
of  immediate  interest  are  the  initial  shock  velocity  Ug  and  the  flow  velocity 
u?  behind  the  shock.  The  flow  velocity  Up  is  also  the  velocity  of  the  inter¬ 
face  or  contact  surface. 

Tine  pressure  and  particle  velocity  are  the  same  on  each  side  of  the  con¬ 
tact  surface.  Assuming  the  driver  gas  expands  isentropically  at  a  constant  7, 
the  speed  of  sound  of  the  driver  gas  at  the  contact  surface  can  be  found  from 


a. 


G 


3> 

V 


V1 


V 


/Pp\  2 7. 

=  hr)  =  (P~P- 


21  14 


% 


(4) 


The  initial  shock  reflects  from  the  nozzle  bringing  the  flow  essentially 
to  rest.  The  velocity  of  the  reflected  shock,  UgR,  and  the  values  of  the 
thermodynamic  properties  behind  the  reflected  shock  are  also  obtained  from 
Reference  9- 
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E  <  ”3  +  P25 
32  -  aa  +  P25 

'I*ue  use  of  a  reflected  rarefaction  vave  is  given  by 


P75  *  f(P75  "  P25^  V  a,p^  +  pj  -  g  -  1  -  0 


r  a  +  p  -| 

vhere  f’  [  J^srA-] 


(10b) 


r  Pp  is 

g « (1  -  p25) 


(10b) 


\5>  3l^§5 

2  r25 


(10c) 


Relations  (8c)  and  (10c)  determine  the  type  of  reflected  vave.  Equations 
(7)  or  (9)  can  be  used  to  calculate  the  pressure  ratio  across  the  reflected 
vave.  'ibis  pressure  ratio  can  be  used  to  calculate  the  flov  velocity  after 
the  vave  from  the  folloving  equations'^ 

U7  1  f  6 

a  =  M5  +  [_^75  ”  *’or  a  rarefac'tion  vave  (ll) 


7  8  M,  + 

a5  5 


P79  ij 

7  Vp(a  p?5+i) 


for  a  shock  vave 


In  using  equations  (7)  through  (12)  the  assumption  is  made  nhat  the  value 
of  7  for  each  gas  does  not  change  across  the  reflected  vave. 

A  special  case  of  this  interaction  betveen  the  shock  vave  and  the  contact 
surface  occurs  when  the  following  relation  is  satisfied 

K  -  g3  *  ?25 
*15  “  +  P 


(13) 


>1 


CONVAIN-SAN  DISOO  CONVAIR  DIVISION  QQ 
OCNKRAL  DYNAMICS  CORPORATION  'Aj 


PAGE  9 

THEO  AERO  NOTE  33 
DATS  6/10/60 


This  is  the  tailored  interface  condition  in  which  there  is  no  reflected  wave. 
Since  there  is  no  reflected  wave  from  the  contact  surface  the  following  relar- 
ticas  must  also  be  satisfied.1 

P3  *  P5 

(HO 

uq  =  u5 


so  therefore 


P0  P5 

P  53  P 

3  2 


u5  -  “a  -  u3  -  u3 


the  change  in  velocity  across  a  shock  is  given  by 


Au 

-  S z 

a 


2 _ /  2  -  1 

TTV-ijJ I  1  *  ^5 


Combining  equations  (15)  and  (l6)  gives 

72(72  -  1) (l  +  a2P)  73(73  -  1)  (1  +  a  D 


which  must  be  satisfied  if  there  is  to  be  no  reflected  wave. 

The  driver  gas  expands  to  the  conditions  in  region  3  through  two  unsteady 

clx 

expansions  and  one  steady  expansion.  The  slope  of  an  expansion  wave,  ^  ,  in 
the  x,t  plane  is  given  by 
cLx 

—  «  u  -  a  for  a  left  running  wave  (l3) 

cLx 

—  *  u  +  a  for  a  right  running  wave  (19) 

In  order  to  construct  the  wave  diagram  it  is  therefore  necessary  to  know  the 
velocity,  u,  and  the  speed  of  sound,  a,  at  each  point  of  interest.  For  a  non 
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a„  +  u  = 


a,  +  u. 


7  -  1  i  i  7  -  1  j  j 
For  a  steady  expansion 


2  2.2  2  .  2,2 

- r-  a,  •!•  u.  = - T  a.  +  u. 

7-1  i  i  7-1  J  j 


and  since  all  expansions  are  adiabatic 
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(20) 


(21) 


(22) 


Equations  (18)  through  (22)  can  be  used  to  calculate  the  flow  velocity 
and  speed  of  sound  at  any  point  in  the  initial  expansion  wave  flow  field. 

The  last  wave  of  interest  in  constructing  a  wave  diagram  to  determine 
steady  flow  duration  is  the  rarefaction  wave  reflected  from  the  end  of  the 
driver.  The  propagation  of  this  wave  is  calculated  using  the  flow  properties' 
determined  from  the  initial  expansion  wave  analysis  described  above.  This 
reflected  rarefaction,  wave  i~  a  right  running  wave  so  the  velocity  of  propaga¬ 
tion  is  given  by  equation  (19)  • 


RESULTS 

The  above  method  of  constructing  wave  diagrams  has  been  applied  to  the 
Convair  Hypersonic  Chock  Tunnel  for  two  sets  of  initial  conditions  using  the 
theoretical  initial  shock  wave  data  of  Figure  k.  One  of  these  initial  condi¬ 
tions  is  for  a  relatively  low  diaphragm  pressure  ratio,  »  5kQ,  and  the 
other  is  for  a  tailored  interface  condition  of  higher  pressure  ratio.  The 
exact  initial  conditions  and  the  resulting  wave  diagrams  are  presented  in 
Figures  5  and  6.  These  wave  diagrams  have  been  computed  for  both  the  present 
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driver,  length  5*4  ft.,  and  for  a  driver  extended  to  13.1  ft.  00  that  tvo 
reflected  rarefaction  vreive3  (ire  slicm  in  each  figure.  For  clarity  the  effect 
of  the  rarefaction  i/ave  reflected  from  the  end  of  the  present  driver  upon  the 
other  waves  in  the  shock  tube  has  not  been  shown.  It  is  readily  seen  in  Fig¬ 
ures  -j  and  6  that  the  reflected  rarefaction  wave  in  the  present  driver  will 
attenuate  the  reservoir  pressure  and  therefore  lead  to  a  time  gradient  of 
dynamic  pressure  at  the  model. 

Figure  7  shews  steady  flow  duration  as  a  function  of  3hock  Kach  number 
with  cut  off  lines  for  three  driver  lengths.  It  will  be  seen  that  the  present 
driver  length  reduces  the  steady  flow  duration  for  Hg  -  9*7*  Tin  effect  of 
this  in  reducing  the  operating  range  of  the  tunnel  is  shown  in  Figure  8.  It 
appears  that  steady  simulation  cf  conditions  below  15,GOO  ft ./sec.  would  be 
eliminated. 

La  Figure  7  the  tailored  interface  condition  is  show?  as  being  achievable 
over  only  a  very  narrow  band  of  shock  Mach  numbers.  In  practice,  however,  the 
shock  waves  and  expansion  waves  reflected  from  the  contact  surface  for  a  band 
of  shock  Macs  numbers  near  this  condition  should  net  be  sufficiently  strong 
to  off  .ct  the  flew  over  the  model  significantly.  Thei’efore,  the  tailored 
interface  concept  may  be  effective  if  good  control  of  shock  velocity  is  attain¬ 
ed. 

If  an  3  ft.  extension  is  made  to  the  driver,  the  nuxinmm  tailored  inter¬ 
face  flew  duration  may  jo  obtained  as  well  as  virtually  eliminating  the  limits 
c>n  the  reflected  shock  case . 

A  decision  on  the  length  of  extension  to  be  selected  should  be  based  on 

* 

flew  d'lrution  measurements  in  the  facility  since  at  least  two  uncertainties 
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apply  to  these  c.alculations.  The  speed  of  sound  in  the  driver  gas  lias  been 
obtained  by  thermodynamic  calculations  of  the  combustion  process  ana  may  be 
in  considerable  error.  Also  no  allowance  has  been  mode  for  shock  attenuation 
due  to  viscous  effects  in  the  driven  section. 
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